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ABSTRACT 

We present an analysis of regular timing observations of the high-magnetic-field 
Rotating Radio Transient (RRAT) J1819— 1458 obtained using the 64-m Parkes and 
76-m Lovell radio telescopes over the past five years. During this time, the RRAT 
has suffered two significant glitches with fractional frequency changes of 0.6 x 10~^ 
and 0.1 X 10"^. Glitches of this magnitude are a phenomenon displayed by both 
radio pulsars and magnetars. However, the behaviour of J1819— 1458 following these 
glitches is quite different to that which follows glitches in other neutron stars, since 
the glitch activity resulted in a significant long-term net decrease in the slow-down 
rate. If such glitches occur every 30 years, the spin-down rate, and by inference the 
magnetic dipole moment, will drop to zero on a timescale of a few thousand years. 
There are also significant increases in the rate of pulse detection and in the radio pulse 
energy immediately following the glitches. 
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1 INTRODUCTION 

Rotating Radio Transients (RRATs) are sporadic radio 
sources from which we detect bursts of emission of typi- 
cally a few milliseconds duration and which occur at inter- 
vals of between a few minutes and a few hours. Eleven of 
these elusive objects were discovered in a single-pulse anal- 
ysis of the Parkes Multibeam Pulsar Survey (McLaughlin et 
al. 2006), but they cannot be detected through their peri- 
odicity using standard Fourier or periodogram techniques. 
However, careful studies of the times of arrival of the bursts 
reveal underlying periodicities which indicate that they are 
rotating neutron stars. Since their discovery, a campaign of 
timing observations has been implemented in order to moni- 
tor the rotational history of these objects (Mclaughlin et al. 
2009). Rotating once every 4.26 seconds, J1819— 1458 is the 
most prolific bursting RRAT. At 1.4 GHz it exhibits ~ 3 ms 
bursts with peak flux density as strong as ~ 10 Jy with a 
burst rate of ~ 20 — 30 hr~^. The high burst rate enabled a 
quick determination of a timing solution for this source and 
hence an accurate position. The true nature of this source as 
a neutron star was further illustrated by X-ray observations. 
These revealed a thermal X-ray spectrum which is what is 



expected from a cooling neutron star (Reynolds et al. 2006; 
Gaensler et al. 2007; McLaughlin et al. 2007). 

The discovery of RRATs raised interesting questions 
about the number of neutron stars in the Galaxy. Initial 
estimates indicate that RRATs may outnumber the nor- 
mal radio pulsar population in the Galaxy by a factor 
of 4 (McLaughlin et al. 2006). It has been argued that 
the observed Galactic supernova rate seems now to be in- 
sufficient to account for the whole population of neutron 
stars (Keane & Kramer 2008), which include radio pul- 
sars, millisecond pulsars, magnetars, RRATs, X-ray de- 
tected but radio quiet Isolated Neutron Stars (INSs) and 
the Central Compact Objects (CCOs), some of which are 
thought to be neutron stars. Evolution between the var- 
ious populations of neutron stars could go some way to- 
wards resolving this problem. Indeed, there is evidence of 
an evolutionary trend connecting young pulsars with mag- 
netars (Lyne 2004). Among the observed properties which 
suggest that the two populations may be linked are very 
low braking indices ( |Lyne et al. 1996 1, magnetar-like burst 
activity in one young pulsar l|Kumar fc Safi-Harb 200^ 
IGavriil et al. 2008]| , low luminosity during quiescence in var- 
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ious magnetars, resembling the X-ray emission of normal 
pulsars ([Archibald et al. 2008|) , and pulsed magnetar radio 
emission (Camilo et al. 2006; Camilo et al. 2007). 

Long-term timing of all these objects provides essen- 
tial information needed to understand any evolutionary re- 
lationships between the different neutron star populations. 
Among the different trends and features observed in long- 
term data we find glitches, which are discontinuous events 
resulting in spin-up of the neutron star, and thought to be 
caused by the transfer of angular momentum from the ro- 
tating superfluid in the stellar interior to the outer crust. 
Glitches sometimes dominate the long-term spin evolution 
of pulsars (Lyne et al. 1996) and have been observed in nu- 
merous young pulsars and several magnetars, appearing to 
be a normal phenomenon among rotating neutron stars. 

We describe the observations in section 2, before dis- 
cussing the unusual glitch activity observed in RRAT 
J1819— 1458 in section 3. In section 4 we discuss the impli- 
cations of these results as they address to the relationship 
between RRATs and other neutron stars. 



2 OBSERVATIONS 

The observations reported here start with the discovery ob- 
servation in August 1998, followed by an intensive 5-year 
series of recordings starting five years later. The measure- 
ments were carried out using the 64-m Parkes Telescope in 
Australia and the 76-m Lovell Telescope at Jodrell Bank in 
the United Kingdom at a frequency of 1.4 GHz. At Parkes, a 
dual-channel cryogenic receiver was used to receive orthog- 
onal linear polarisations. Each channel was processed in a 
512 X 0.5-MHz filterbank, added in polarisation pairs to pro- 
vide the total intensity, 1-bit sampled and recorded on tape. 
Observations at Parkes were typically 30 min in duration. 
At Jodrell Bank, a dual-channel cryogenic receiver was sen- 
sitive to the two hands of circular polarisation which were 
each processed in a 64 x 1-MHz filterbank before being added 
to produce total intensity, 1-bit sampled and stored on tape. 
Observations at Jodrell Bank were 60 min in duration. The 
longer observations at Jodrell Bank were due to the lower 
sensitivity as a result of the narrower available bandwidth. 
We note also that the presence of broadband impulsive ra- 
dio frequency interference hinders the identification of pulses 
from RRATs at both observatories. However, the quality of 
the observational data was markedly improved by the use of 
the zero-DM filtering algorithm (Eatough et al. 2009). 

Off-line, the data were dedispersed at the nominal dis- 
persion measure (DM) of the RRAT as well as at a num- 
ber of adjacent DM values to aid RFI discrimination. The 
resulting time sequences were then searched for significant 
single bursts of radiation, whose times of arrival (TOAs) 
were recorded. We can compare the observed TOAs with a 
model of the neutron star's rotational and astrometric pa- 
rameters in a manner similar to standard pulsar timing (see 
Lorimer & Kramer (2005) or Lyne & Graham-Smith (2006)). 
The only difference is that standard pulsar methods obtain 
TOAs from integrated pulse profiles of as many as 10^ - 10" 
single pulses whereas for RRATs individual pulses have a 
sufficiently high S/N that we can obtain TOAs from indi- 
vidual pulses (and indeed are required to do so because of 
the sporadic emission). The differences between the observed 
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Figure 1. a) The timing residuals from the TOAs of the indi- 
vidual pulses from J1819— 1458, relative to a simple slow-down 
model, showing that they are located in 3 clearly identifiable 
bands which are offset in pulse longitude by ±45 ms. b) the same 
residuals, with the TOAs in the top band and the bottom band 
decreased and increased by 45 ms respectively. The rms of the 
residuals decreases from 21.2 ms to 9.1 ms. 

arrival times and the predicted arrival times are known as 
the pulsar's "timing residuals" and the model parameters 
(such as period, period derivative and sky position) are de- 
termined by adjusting their values to minimise the rms value 
of the residuals. 



3 RESULTS 
Tri-Modal Residuals 

In Fig. [1] we plot the timing residuals for all the bursts 
received from J1819— 1458 over the most recent 200 days 
relative to a simple slow-down model which includes only the 
spin frequency and its first derivative. The upper plot shows 
that the bursts from J1819— 1458 arrive over a longitude 
range spanning approximately 120 ms and grouped within 
three longitude regions separated by about 45 ms. 

Fortunately, within a single 30-min or 60-min observa- 
tion there is a sufficiently large number of bursts to allow 
each burst to be uniquely identified with one of the three 
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Figure 2. Top: Distribution of the timing residuals of bursts 
from RRAT J1819-1458 where the three bands seen in Fig. [TJi 
are clear. Bottom: A scatter plot of pulse energy versus timing 
residual. 



bands. Those TOAs identified with the upper and lower 
bands are either decreased by 45 ms or increased by 45 ms 
respectively and the residuals recalculated. The lower dia- 
gram in Fig. [1] shows the residuals resulting from this re- 
calculation. As a result of this procedure the rms value of 
the residuals is reduced by a factor of ^ 2.5 and the uncer- 
tainties in the fitted parameters are similarly reduced. The 
timing analyses presented in this paper are based upon all 
the TOAs modified in this fashion. 

Most of the bursts (~ 60%) come from the middle band, 
where also the brightest bursts are found. The top panel of 
Fig.[2]shows a collapsed histogram of the timing residuals of 
the top panel of Fig. [T] This is essentially a probability dis- 
tribution in rotational phase for the bursts. The three-band 
structure is consistent with the three-component profile re- 
ported by Karastergiou et al. (2009). In Fig. O we show a 
sequence of pulses detected during a single long Parkes ob- 
servation, along with the composite profile formed by sum- 
ming all of these individual pulses, which show anywhere 
from one to four components. Such composite profiles vary 
from day to day due to the large degree of pulse-to-pulse 
variation. 

Esamdin et al. (2008) reported the presence of only two 
preferred longitude ranges, using the 25-m radio telescope 
at Urumqi. However, because of the smaller size of that tele- 
scope, they were only sensitive to the most intense bursts. 
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Figure 3. All 165 single pulses detected during an eight-hr obser- 
vation of J1819— 1458 with Parkes at 1.4 GHz (bottom) along with 
the composite profile formed by summing all 165 pulses (top). The 
time intervals between the pulses vary. 

those having fiux densities > 3.4 Jy, resulting in the observa- 
tion of only 162 bursts in 94 hours of observation, or a burst 
rate of 1.7 hr~^. Most of the bursts came from the middle 
band, with 18 from the early band and just 2 from the later 
band. The higher sensitivity of the larger telescopes used in 
the present experiment allowed us to detect > 500 bursts 
from RRAT J1819— 1458 in 27 hours of observation at rates 
of ~ 15 — 20 hr~^ at Jodrell Bank, and somewhat higher at 
Parkes. The bottom panel of Fig. [2] shows the intensity dis- 
tributions of bursts in the three bands. In determining pulse 
intensities, we used the sky temperature model of Haslam et 
al. (1982) extrapolated from 430-MHz to 1400-MHz assum- 
ing a spectral index of —2.6 (Lawson et al. 1987). Taking 
into account the higher flux density limit of the Urumqi ob- 
servations and the unknown amplitude distribution within 
the bands at high fiux density, the two sets of results appear 
consistent. 

The three-component profile is consistent with the tri- 
modal timing residuals, and can be reconciled with a patchy 
emission beam with core and conal components (e.g. Lyne 
& Manchester 1988). We expect that many normal pulsars 
with profiles with multiple components would exhibit similar 
bi- or tri-modal residuals if timed through their single pulses. 

Glitches 

Fig. |4^ shows the evolution of the rotational frequency over 
a ~10-yr interval. In the upper diagram, the steady slow- 
down is interrupted at around MJD 53900 by a sudden 
spin-up. The nature of this spin-up can be seen more clearly 
in Fig. |4]d, which displays the frequency residuals relative 
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Table 1. The observed and derived rotational parameters of the 
RRAT J1819-1458. 
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Figure 4. The frequency evolution of the RRAT J1819-1458 
over a 10-year period, a) shows the secular slowdown in rotation 
rate of the neutron star, interrupted by a major glitch which is 
seen as a clear discontinuity at around MJD 53900. b) shows the 
frequency residuals relative to a simple slow-down model fitted 
to data between MJD 51000 and 53900 and reveals the presence 
of a second, smaller glitch, about 200 days later, c) presents the 
variation in the magnitude of frequency derivative \v\, showing a 
significant decrease in the rate of slowdown following the glitches, 
c) shows the timing residuals trcs relative to the rotational model 
given in Table 1. 



Timing parameters 


Right Ascension o 


18''19™34M73 


Declination 5 


-14°58'03".57 


Frequency v (Hz) 


0.23456756350(2) 


Frequency derivative i> (s~^) 


-31.647(1) X 10-15 


Timing Epoch (MJD) 


54451.0 


Dispersion measure DM (cm~'^pc) 


196 


Timing data span (MJD) 


51031 - 54938 


RMS timing residual a (msec) 


10.2 


Glitch 1 Parameters 


Epoch (MJD) 


53924.79(15) 


Incremental Av (Hz) 


0.1380(6) xl0-« 


Incremental Ai> (s~^) 


0.789(6) X 10-15 


Decay &v (Hz) 


0.0260(8) X 10-6 


Decay timescale r (days) 


167(6) 


Glitch 2 Parameters 


Epoch (MJD) 


54168.6(8) 


Incremental Av (Hz) 


0.0226(3) xlO-'' 


Derived Parameters 


Characteristic Age (kyr) 


120 


Surface Magnetic Field (G) 


50 X 10^2 



we used the position of J1819— 1458 derived from Chandra 
X-ray observations of the source (Rea et al. 2009), which 
has errors of 0.2 arcsec in both coordinates. The value of 
dispersion measure used was that given by McLaughlin et 
al. (2006). 

Glitches usually display fractional increases in rota- 
tional frequency ^vjv which range between 10"^ and 10"^ 
(e.g. Lyne, Shemar & Graham-Smith 2000). The two glitches 
in J1819-1458 have ^^ujv = 0.6 x IQ-*^ and 0.1 x IQ-^ 
While these are not the largest glitches known (c.f. Hobbs 
et al. 2002) they are typical in fractional size to those from 
many young pulsars (e.g. Zou et al. 2008) and magnetars 
( |Mereghetti 2008[ ). 



to a simple slow-down model fitted to the TOAs obtained 
prior to MJD 53925. The spin-up is revealed as being due to 
two glitches, quite close to one another, occurring at around 
MJD 53926 and MJD 54167, followed by a period of steadily 
increasing frequency relative to the model. The latter corre- 
sponds to a reduced magnitude of frequency first derivative, 
which can be seen in Fig.|4}; following a transient short-term 
increase. 

In order to quantify these changes, the TEMPCQ pulsar 
timing package was used to model the slowdown of the pulsar 
and the two glitches. It was found necessary to include steps 
in both V and v and a decaying exponential for the first glitch 
but only a step in v for the much smaller second glitch. These 
parameters are given in Table 1. For this timing analysis, 

1 http:/ /pulsar. princeton.edu/tempo 



Glitch-associated emission activity 

In Fig. [S] we show the rate of pulse detection (top) and 
the peak pulse energy (bottom) for each observation as a 
function of date for J1819— 1458. For both quantities, there 
is a large variation, with maximum values in the observation 
(on MJD 53960) which immediately follows the first, and 
largest, glitch. On this date, the mean pulse detection rate is 
68±12 pulses/hour, 2.8 times the mean and significant at the 
3.5(T level. The corresponding value of the peak pulse energy 
is 58.7 mjy.s, which is 3.7 times the mean and significant at 
the 4.7a level. We note, however, that the high burst rate 
on MJD 54619 is not associated with any obvious timing 
abnormality. 
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Figure 5. The variation of rate of pulse detection (top) and 
the peak pulse energy (bottom) with date for J1819— 1458. The 
epochs of the two glitches are marked with dashed lines. As a con- 
trol experiment, we find no dependence of the pulse detection rate 
on varying sensitivity due to the zenith angle of the observations. 



4 DISCUSSION 

While J1819— 1458 exhibits many of the properties of nor- 
mal, relatively young radio pulsars, namely rapid spin-down 
and glitches, the behaviour following the glitch activity is 
quite different. The main effect is a long-term reduction in 
the magnitude of the slow-down rate. This is unlike all ob- 
served glitches in radio pulsars and magnetars, which always 
lead to a long-term increase in the rate of slow-down. Fig. [6] 
illustrates this by showing the behaviour of slow-down rate 
for J1819— 1458 and a selection of normal young pulsars for a 
few hundred days around typical glitches. The behaviour of 
the Crab pulsar (B0531+21) shows the greatest similarity 
to that of J1819— 1458, each displaying a relatively short- 
term transient before reaching a new asymptotic slow-down 
rate, although the long-term net increments in slow-down 
rate clearly have different signs. 

Fig.[7]shows the P-P diagram for normal pulsars, mag- 
netars and other RRATs QMcLaughlin et al. 2009| |. Objects 
with constant dipolar magnetic field move towards the lower 
right-hand corner of the diagram with a slope of —1. Young 
normal pulsars are observed to move with a slope of between 
— 1.0 and -fO.5. On this diagram, glitches from normal pul- 
sars such as those shown in Fig. |S] result in a net increase in 
slow-down rate and an upwards step in the P-P diagram. 
On the other hand, RRAT J1819-1458 stepped vertically 
downwards, towards smaller values of P. If this particular 
post-glitch behaviour is typical, then the long-term effect of 
any glitches would be a secular movement towards the bot- 
tom of the P-P diagram. If such glitches were to occur every 
30 years, then the slowdown rate would decay to zero on a 
timescale of only a few thousand years. Only larger time- 
span observations will unveil the actual path of RRATs on 
the P-P diagram. If the trend continues it could indicate 
that the RRAT started off in the region of the diagram pop- 
ulated by the magnetars. 

Glitches are understood to be caused by the communi- 
cation between a superfluid component of the neutron star 
interior and the crust. In the standard model, a neutron star 
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Figure 6. The changes in the magnitude of the first derivative 
of the rotational frequency |!>| near the glitch from the RRAT 
J1819— 1458 and near typical glitches in 11 pulsars. A second 
frequency derivative has been fitted to each pre-glitch data set 
and subtracted from the whole set. The scale of each diagram has 
been adjusted so that the step change in value at the glitch is the 
same. Without exception, the final values of for the pulsars 
are greater than before the glitches. In the case of J1819— 1458, 
there is a net decrease in magnitude. 



has a superfluid that is originally spinning faster than the 
crust (see, e.g. Anderson & Itoh 1975, Alpar & Pines 1993). 
This is because the magnetic dipole acts on the crust and the 
coupled core, with the crustal superfluid rotating indepen- 
dently. All neutron stars, regardless of their observational 
manifestation, are expected to show this same structure. 
The angular momentum resides on vortex lines, which are 
pinned to nuclei in the crust. These lines suffer strong forces 
due to the mismatch between the superfluid and crustal ve- 
locities. A glitch occurs when there is a sudden unpinning of 
the vortex lines, and subsequent transfer of angular momen- 
tum from the superfluid to the crust. In a magnetar, glitches 
could occur due to the high internal magnetic fleld that can 
deform or crack the crust ( Thompson fc Duncan 1996^ . 

It could be that the strange glitch behavior of 
J1819— 1458 is due to a different mechanism for the un- 
pinning, namely the deformation or cracking of the crust 
due to the high magnetic field, as opposed to the unpin- 
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Figure 7. The P-P diagram for normal young pulsars and mag- 
netars, as well as for the RRAT J1819— 1458 (filled diamond) and 
other RRATs (open diamonds). 



ning due to the angular velocity lag as suggested for normal 
pulsar glitches. While a magnetar glitch which resulted in a 
decrease in the spin-down rate has not yet been observed, 
magnetar glitches seem to show different characteristics to 
radio pulsar glitches. Gavriil et al. (2009) showed that AXP 
4U 0142+61 suffered a glitch that resulted in a net spin- 
down. They interpret this as being due to some regions of 
the superfluid originally spinning slower than the crust, as 
opposed to faster in normal rotation-powered pulsars. 

Magnetar glitches are occasionally associated with ra- 
diative events and pulse profile and spectral changes (see, 
e.g. Dib et al. 2008), with no obvious relationship between 
the size of the glitch and the extent of these changes. The 
increase in both burst rate and peak pulse energy immedi- 
ately following the first glitch in J1819— 1458 is suggestive of 
an association. Both these quantities seem to be sustained 
at a raised level for 100-200 days. Although there are some 
days with nearly as high a burst rate that are not associated 
with any timing irregularity, the fact that both burst rate 
and pulse energy are associated with the glitch suggests that 
this is not just a statistical fluke. More large glitches must be 
detected for such correlations to be tested robustly. In ad- 
dition. X-ray observations immediately after the next large 
glitch would be useful to search for magnetar-like bursts or 
correlated pulse proflle or spectral changes. 



5 CONCLUSION 

In this paper, we have demonstrated the difficulty of tim- 
ing pulsars with multiple components through their single 
pulses, and have presented a simple algorithm for mitigating 
these effects. 

We have presented a timing analysis of the RRAT 
J1819-1458 over the past 5yr and report the detection of 
two significant glitches. The magnitudes of these glitches 
are similar to those measured for radio pulsars and mag- 
netars, but the behaviour after the glitch is quite different. 
The glitches have resulted in a long-term net decrease in the 
slow-down rate, corresponding to a downward vertical move- 
ment on the P — P diagram. This jump in frequency deriva- 
tive is ~ 0.75(5) X 10~^^ Hz s~^, corresponding to a ~1% 
decrease of the surface magnetic field of the neutron star. 
If such glitches are representative behaviour of this RRAT 
and were to occur every 30 years, the spin-down rate, and 
by inference the magnetic dipole moment, will drop to zero 
on a timescale of a few thousand years. This implies that 
J1819— 1458 may have begun its life in the magnetar region 
of the diagram. 

We have also briefly discussed the evolution of the burst 
rate and peak flux density of pulses from J1819— 1458 with 
time, and showed that both peaked immediately after the 
large glitch. This is tantalizing evidence for a similarity with 
magnetar glitches. 
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